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ABSTRACT 

Using the Caltech Submillimeter Observatory 10.4-meter telescope, we per- 
formed sensitive mapping observations of ^^CO J = 7 — 6 emission at 807 GHz 
towards Orion IRc2. The image has an angular resolution of 10^', which is the 
highest angular resolution data toward the Orion Hot Core published for this 
transition. In addition, thanks to the on-the-fly mapping technique, the fidelity 
of the new image is rather high, particularly in comparison to previous images. 
We have succeeded in mapping the northwest-southeast high-velocity molecular 
outflow, whose terminal velocity is shifted by ^ 70 — 85 km s~^ with respect to 
the systemic velocity of the cloud. This yields an extremely short dynamical time 
scale of ^900 years. The estimated outflow mass loss rate shows an extraordi- 
narily high value, on the order of 10~^ Mq yr~^. Assuming that the outflow is 
driven by Orion IRc2, our result agrees with the picture so far obtained for a 20 
Mq (proto)star in the process of formation. 

Subject headings: ISM: jets and outflows — ISM: molecules — stars: early-type 
— stars: individual (Orion Hot Core) — submillimeter 
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1. Introduction 

Despite its astrophysical importance, the high-mass (M* ^ 8 Mq) star formation pro- 
cess remains poorly understood both observationaUy and theoreticaUy. High-velocity (HV) 
molecular outflows are the most prominent phenomena seen, not only in high-mass but also 
in low-mass star-forming regions. Molecular outflows from low-mass stars are thought to be 
momentum- driven by highly collimated stellar winds, and are closely related to the accretion 
process onto the central star (e.g., Arce et al. 2007, and references therein). Although such 
a paradigm cannot simply be scaled up to the high-mass regime, one may obtain important 
information about the formation mechanism of early- type stars by observing their outflows. 
In this context, the Orion Hot Core (HC) is the best site to study massive star formation 
process, as its proximity allows very detailed studies. 

Ground-based mapping observations of CO lines at s u bmill imeter (submm) wavelengths 



tow ards the region were p ioneered by, e.g., I Wilson et al.l ( 120011 ) in the J = 7 — 6 transition 



and iMarrone et al.l (120041 ) in J = 9 — 8. These authors focused on studying the properties 
of the intense low-velocity gas around the velocity of the natal molecular cloud, which arises 
predominantly from a photo-dissociation region (PDR). Observations of low- J lines, such 
as J = 2 — 1, allow observers to attain high angular resolution imag e s down to by 



utilizing interferometers, as demonstrated by IChernin fc Wrightl (119961 ). iBeuther fc Nissen 



( 120081 ). and others. However, such low- J line observations of molecular outflows usually make 
it diflicult to distinguish low- velocity outflowing gas from the quiescent ambient gas. The 
J = 7 level is approximately 156 K above the ground-state level, and thus more sensitive to 
the presence of hot (T > 100 K) gas. The Einstein A coefl&cient for the J = 7 — 6 transition 
is 50 times greater than for the J = 2 — 1 line, providing a high contrast between the warm 
and cold gas. Such characteristics greatly help us to distinguish HV outflowing gas from 
the surrounding material. To characterize the physical properties of the outflowing gas in 
the Orion HC, we performed observations of CO J = 7 — 6 using the Caltech Submillimeter 
Observatory (CSOj^ 10.4-meter telescope, making full use of its high spatial and high spectral 
resolution. 



-^Caltech Submillimeter Observatory is operated by the California Institute of Technology under the grant 
from the US National Science Foundation (AST 08-38261). 
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2. Observations and Data Reduction 

We carried out on-the-fly (OTF) mapping observations of the CO J =7^6 line (rest fre- 
quency, z^rest = 806651.720 MHz) with the CSO on 2009 January 21 and 22. The observations 
took 1-2 hrs each day. Prior to the submm observations, we performed optical pointing ob- 
servations on 2009 January 20. Therefore, we believe that the overall pointing accuracy was 



better than ?>" . We used the 850 GHz receiver (IKooi et al.ll2000l ) with the Fast Fourier Trans 



form Spectrometer (FFTS) as a back-end, yielding an effective velocity resolution (A'^res) of 
0.0454 km s~^ with the 1 GHz bandwidth mode. During the observations, the atmospheric 
zenith optical depth at 225 GHz ranged between 0.035 and 0.055, and the single side-band 
(SSB) system temperature (Tgyg) stayed between ^2500 K and 8000 K, depending on the 
optical depth and the airmass. We wish to point out th at our observa t ions a re almost twice 



as sensitive as the previous ^^CO (7-6) observations by I Wilson et al.l (120011 ). whose typical 
Tgys was ^ 14, 000 K. The OTF mapping was carried out over an area of ^ 170^' centered 
on Orion IRc2 (R.A. = 5^35^14.5^ Dec. = -5° 22^ 30^.^4 in B1950) by scanning once each 
along the R.A. and Dec. directions. Subsequently, we decided to concentrate on mapping 
the central 85^' region, which we further scanned twice along the R.A. and once along the 
Dec. directions. The OTF mapping scans were gridded with a pixel size scale of 5^.'0 by 
employing a position-switching method. The telescope pointing was checked by observing 
the continuum emission towards Venus and Saturn every hour. At the line frequency, the 
beam size (^hpbw) is estimated to be 10^' and the main-beam efficiency {vj^h) was 0.34 from 
our measurements towards Saturn. All the spectra were calibrated by the standard chopper 
wheel method, and were converted to main-beam brightness temperature (T^b) scale by di- 
viding by rymb- The uncertainty in the intensity calibration is estimated to be ^20%. After 
calibrating all the spectra, we made two 3D data cubes from two groups of the spectra taken 
by scanning along R.A. and Dec. directions. These two cubes were processed with the code 
"Basket- Weave" , implemented in the NOSTAR package (Sawada et al. 2007) to remove the 
scanning effect using the method developed by Emerson & Graeve (1988). 



3. Results and Discussion 



Figure [T] shows our CO (7-6) spectrum towards Orion IRc2 in the T^b scale. The 
spectral profile is rather broad, and shows prominent HV wing emission. The HV wing 
is seen both blue- and redshifted sides up to ~80 km s~^ away from the systemic velocity 
(^ys) of the cloud, yLSR= 9 km s ~^. The peak T^h is 202 K at I^lsr^ 1-9 km s~-^, which is 
consistent with that measured by IComito et al.l (120051 ) , although the pointing centers differ 
between our maps and theirs by 5'.'0. The CO spectrum shows a dip at l^yg, suggesting that 
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the line is optically thick around l^yg. The two isolated emission lines can be seen at ^lsr"^ 
+110 km s~^ and ^ —115 km s~^; they are identified as lin es of NS (frpst, =80 6350 MHz) 
and ^^CHaOH (806973 MHz), on the basis of the line survey dComito et al.lbood ) at the 350 
jam band. 

Figure [2] shows an o verlay of the low - velocity bulk emission of the CO (7^6) and the 350 
/im continuum emission (ILis et al.lll998l ). The (7-6) line shows intense emission surrounding 
IRc2, and spreads along the north-south direction if we consider the weak extended compo- 
nent. The CS Q beam is to o large to resolve the spatial separation of '^0'.'8 between IRc2 
and Source I (IGezaril 1 19921 ) which has been proposed as a powering source for the Orion 
outfiow (Beuther & Nissen 2008, and references therein). The north-south elongation of the 
weak emission corresponds to the Orion Ridge; the northern tip of the Orion South core can 
be recognized at the bottom (i.e., south) of our map. The widespread quiesc ent CO (7-6) 



emiss ion is thought to be associated with the PDR heated by Q^C Orionis (j Wilson et al. 
200 ll ). The bulk emission shows a condensation with a diameter of ^ 80^', measured at the 
50% level contour, corresponding to 0.065 pc at = 450 pc. It should be noted that the CO 
(7^6) line and the 350 /xm continuum emission do not show similar spatial structure; the 
former has a roundish shape, whereas the latter is elongated along the north-south direction. 
Another result deduced from Figure [2] is that the bulk emission does not show a single peak 
at the position of the 350 /xm continuum peak. Instead, it shows two local maxima; one of 
them lies to the west of the BN object and the other is ^ 10^' east of IRc2. Positions of the 
two local maxima are roughly consis tent with the peak p ositions of the blue- and redshifted 
CO (7-6) wing emission reported in I Wilson et al.l (120011 ). although their HV wing emission 
maps (see their Figure 5) show a scanning eflFect. 

Figure [3] shows velocity channel maps of the HV emission. To produce the channel maps, 
we convolved the data with a Gaussian beam so that the effective beam size (^efr) becomes 
13^', to obtain a higher si Rnal-to-noise r atio ( S/N). Here, the choice of 13^' is to compare 
with the previous work by I Wilson et al.l (120011 ). who used the HHT 10-meter telescope. The 
channel maps show that the redshifted gas is elongated to the (south)east of the central star, 
whereas the blueshifted gas seems to be distributed westward. 

The distribution of the blue- and redshifted HV wing emission is compared in Figure H] 
by overlaying the integrated intensity maps. Clearly, the HV gas is confined to the innermost 
region of the core, inside the 50% level contour of the low-velocity bulk emission. Here, we 
integrated the HV emission over the velocity ranges shown in Figured! i.e., over the velocity 
range between the boundary velocity (H), which divides the LSR- velocity range of the out- 
fiowing gas from the bulk emission, and the first highest LSR- velocity where the wing drops 
below the 5a level. We refer to such an LSR- velocity as the terminal velocity (Vt). Since it is 
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almost impossible to separate low-velocity outflowing gas from the bulk ambient gas, we ar- 
bitrarily estimated 14 from Figures [1] and [3l The redshifted gas shows an extended structure 
along the east-west direction, whereas the blueshifted emission arises from a compact con- 
densation, rather than an elongated structure. It appears that the distribution of the blue- 
and redshifted gas is symmetric with respect to the peak position of the 350 fim continuum 
source, although the blue- and redshifted gas overlap each other. We pre fer to interpret this 
as rep r esenting the pair of molecular o utflow lobes prev i ously r eported [Chernin fc Wright 
(Il996l ): iRodrfguez-Franco et all fllQQoh : Ide Vicente et all J^OoJ; ISeuther fc NissenI fl2008l y|^ 
This is because their terminal velocities are too high to be interpreted in the context of 
the other scenarios, such as rotation or/and infall. Given its position, we exclude the BN 
object as a candidate source for the outflow. However, due to the limited angular resolution, 
our data do not give furthe r constraints on identifyin g the driving source, as discussed by 
e.g.. Ide Vicente et al.l (120021 ) . iBeuther fc NissenI (120081 ) . iRodriguez et al.l (120091 ). and others. 
Assuming that the position of IRc2 represents the center of the outflow, we made a position- 
velocity (PV) diagram of the CO (7^6) emission along the outflow axis (Figure [5j). The PV 
diagram principally conflrms our flndings described above. Unfortunately, the CSO beam is 
too large to provide more detailed information about, e.g., driving mechanisms for outflow 
and/or outflow history. 

Subsequently, we estimated the outflow lobe mass (Miobe) to calculate kinematical prop- 
erties such as mass-loss rate (Mflow) and momentum rate (Fco). Here, we assumed that the 
HV gas is optically thin and in LTE because neither isotope line, i.e., ^^CO J = 7 — 6, nor 
similar high- J lines, e.g., ^^CO J = 6 — 5, are available with the 10^' beam. We obtained 
M\ohe = (1.3 ±0.3) and (3.5 ±0.9) Mq for the blue and red lobes, respectively, by obtaining a 
mean column density over the lobes of (9.0 ±2.4) x 10^^ and (1.0±0.3) x lO"*^^ cm~^. Here, we 
adopted an excitation tem perature (Tpx) o f 150 K (IComito et al.ll2005l ). and used a ^^C0/H2 
abundance ratio of 10~^ (jPickmanl 1 19781 ). For comparison with previous publications, we 
adopted d = 450 pc, although recent astrometry using the VLBI technique has revised the 
distance to the region (414±7 pc, Menten et al. 2007; 437±19 pc, Hirota et al. 2007). It is 
likely that the outflow masses are underestimated because of the uncertainty in deflning V]j. 
Since the highest outflow velocity is given by Vnow = \Vt — l^ys|, the dynamical time scale (t^) 
is estimated from /lobe/^ow, where the lobe length /lobe is deflned as the maximum extent of 
the lobe measured from the IRc2 position. We assumed that the outflow inclination angle (i) 
is 45°, here deflned as the angle between the outflow axis and the line of sight. The estimated 
^ 900 years suggests that the outflow is extremely young. Hereafter, we describe the mean 
values of the outflow properties, because the derived outflow properties for the blue- and red- 
shifted lobes are comparable within the errors. We can estimate the outflow mass-loss rate 
from Miobe/^d- This is ^ 4 x 10~^ Mq yr~^ for both lobes. Since molecular outflows appear 
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to be momentum-driven (jCabrit fc Bertoutlll992l ). the momentum rate Fco=Mobe^ow/^iobe 



may be taken as an indicator of the outflow strength and hence of the mass and luminosity 
of the young steUar object (YSO) powering it. The outflow lobe has Fco ^ 0.3 Mq km s~^ 
year~^. The derived properties are not aflFected signiflcantly if one takes into account the 
unknown inclination angle which is deflned as the angle between the outflow axis and line 
of sight. In fact, td, ^flow, and Fco are corrected with factors of cotz, tani, and sini/ cos^i, 
respectively (see e.g., Davis et al. 1997). Thus, even if we assume extreme values of z ^ 20° 
and 70°, the corrections for the three quantities above would only be a factor of 3-5, thus 
leaving our estimates unaffected. All the derived outflow parameters are extraordinarily 
large, implying that the powering YSO has a bolometric luminosity of ^ 10^ — 10^ L©, if 
we apply the empirical relationship between the two quantities (see Figure 5 of Richer et al. 
2000 and Figure 4 of Beuther et al. 2002). 

Our estimate of Mflow to be on the order of 10~^ Mq yr~^ strongly suggests that the 
central massive YSO is gaining mass toward its flnal stellar mass through a disk-outflow 
system. We do not exclude the possibility that mass accretion onto the central star (or/and 
disk system) is currently enh anced. A stat i stical study of molecular outflows in 26 high- 
mass star- forming regions by iBeuther et al.l (12002 ) suggests that the ratio between the jet 
mass-loss rate and mass-accretion rate (Mace) be ^ 0.2. If this is the case, the accretion 
rate onto the outflow powering source would be on the order of ^ 10~^Mq yr~^. Although 
the driving source of the Orion outflow has not been flrmly identifled, and has been a 
matter of debate (de Vicente et al. 2002; Beuther and Nissen 2008; Rodriguez et al. 2009), 
it is tempting to_suggest that this accretion rate does not contradict that estimated by 
Nakano et al.l (120001 ). who proposed that the central massive YSO in Orion IRc2 is either 
a protostar with a stella r radius of 30 Rq or a protostar which has a disk with 



10-2 Mq yr- 



Nakano et al.l (120001 ) used a one-zone model to account for the 
observational result s of R^ ^ 300 i?© and an eflFective temperature of T^^^ 3000 - 5500 K 



(IMorino et al.lll998l ): they adopted the latter interpretation because their model could not 



reproduce an Lboi ^ 4 x 10^ (proto)star (Kaufman et al. 1998; Gezari et al. 1998) using 
the low Teff a nd a stellar radius la r ger th an 30 Rq. More stringent calculations recently 
performed by iHosokawa fc Omukail (120091 ) considered the possibility that such a massive 
protostar can be formed if Mace is higher than 4 x 10~^ Mq yr~^. If this is the case, the 
Keplerian speed at the protostellar surface ('^kep = \/ GM^/R^) , which is believed to represent 
the jet velocity (e.g.. Richer et al. 2000; Arce et al. 2007), becomes ^ 220 km s~^ for a 25 



Mm an d 100 Rq star. Here, the stellar mass and radius are taken from iHosokawa fc Omukai 
(120091 ). We believe that our result of Vqow ~ 80 - 90 km s~^ agrees with the picture of 
Orion IRc2 so far obtained, especially when we consider that molecular outflows are made 
of ambient material entrained by a jet/winds from the central star. In conclusion, our 
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sensitive and high-fidelity ^^CO (7-6) hne imaging demonstrates that a massive YSO is 
driving an extraordinarily powerful molecular outflow, suggestive of a large accretion rate. 
Clearly, higher resolution and better sensitivity imaging of the submm transition, as well as 
identification of its driving source, are required to study the spatial structure of the outflow. 
It may retain a history of the outflow events, a clue to understanding the mass accretion 
history. 
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Fig. 1. — Single-dish spectra of the ^^CO 7 — 6 emission toward Orion IRc2 in the main-beam 
brightness temperature (Tmb) scale with a velocity resolution of 0.726 km s~-^. See the y-axis 
labels of the right- and left-hand sides for the overall (filled histogram) and the magnified 
spectra, respectively. The RMS noise level is 0.72 K for 81 seconds integration. The blue- 
and red-colored horizontal bars below the spectra indicate the LSR-velocity ranges used to 
obtain the integrated intensity map shown in Figure [H The systemic velocity of the cloud 
is ^LSR ^ 9 km s~^ which is indicated by the vertical dashed line. 
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Fig. 2. — Overlay of the bulk emission of the -"^^CO (7-6) (green image plus white and green 
contours) on the 350 /xm continuum emission (magenta contours; Lis et al. 1998). The 
contours for the ^^CO (7-6) emission are the 95% level (thick white), 90%, 80%, 70%, 60% 
(thin white), 50% (thick green), and 30% (thin green) level of the peak for the clarity of 
the plot. The CO map is obtained by integrating the emission between 1^lsr= 4.5 km s~^ 
and 14.5 km s~^. The li ght blue triangl e s with source names sh o w th e positions of the 
radio continuum sources Menten fc Reid (j 19951): i Rodriguez et al. (2009)], and the yellow 
stars indicate the positions of the Trapezium members. Note t hat the peak position of the 



infrared source IRc2 is displaced from the radio source I by 0'.'8 ( Gezari 



1992|). The hatched 



ellipse at the bottom left corner indicates the HPBW of the CSO beam (10^^). 
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Fig. 3. — Velocity channel maps of the high-velocity ^^CO (7^6) emission toward the Orion 
IRc2 (the star) with an effective angular resolution of 13^^ (see text). Each channel map is 
averaged over a 5.0 km s~^ bin whose central LSR-velocity is shown on the top of each panel. 
All the contours are 2a intervals, increasing from the 3a level where la is 11 K in Tmb- 
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Fig. 4. — Comparison of the blue- and redshifted -^^CO (7-6) high-velocity (HV) wing emis- 
sion. The green contour, taken from Figure O represents the 50 % level co ntour of the CO 
bulk emission, and the greyscale shows the 350 fim continuum emission (ILis et al.l Il998l ). 
The contours for the HV emission are 2a intervals, increasing from the 4a level where a = 
10.6 and 9.9 K- kms for the blue- and redshifted wings, respectively. The blue- and red- 
shifted emission are integrated over the velocity ranges of —75.5 <l^LSR/km s~^< —30.5, 
and +34.5 <l^LSR/km s~-^< +79.5, respectively (see the horizontal color-coded bars in Fig- 
ured]). The hatched circle at the bottom-right corner indicates the effective beam size (^eff) 
of 13'^ The dashed line indicates the cut axis for the position-velocity diagram shown in 
Figure [51 The other symbol s are the same a s thos e in Figure [21 Notice that the possible 
heating sources discussed by lde Vicente et al.l (120021 ) are located within a region of ^V!5 ra- 
dius, whose spatial-scale is almost comparable with the blue-triangle representing the IRc2 
position. 



-13- 




-60 -40 -20 20_^ 40 60 80 
Vlsr (km s ) 



Fig. 5. — Position-velocity (PV) diagram of the -^^CO (7-6) emission along the outflow axis 
(P.A.= +126°) passing through the position of IRc2 (see Figured]). The white and green 
contours are the same as those in Figure El and the two yellow contours indicate the 20% and 
10% levels with respect to the peak intensity. The vertical dashed-line at V^sr = 9 km s~^ 
shows the Kvs- 



